Mon. Not. R. Astron. Soc. OOP. HICT ('2005') Printed 2 February 2008 (MN I*TeX style file v2.2) 



IGM metal enrichment through dust sputtering 



Simone Bianchi^ and Andrea Ferrara^ 

^ Istituto di Radioastronomia / CNR - Sezione di Firenze, Largo Enrico Fermi 5, 50125 Firenze, Raly 
^ SISS A/International School for Advanced Studies, via Beirut 2-4, 34013 Trieste, Raly 



January 2005 



ABSTRACT 

We study the motion of dust grains into the Intergalactic Medium (IGM) around red- 
shift z = 3, to test the hypothesis that grains can efficiently poUute the gas with metals 
through sputtering. We use the results available in the literature for radiation-driven 
dust ejection from galaxies as initial conditions, and follow the motion onward. Via 
this mechanism, grains are ejected into the IGM with velocities > 100 km s~^; as they 
move supersonically, grains can be efficiently eroded by non-thermal sputtering. How- 
ever, Coulomb and coUisional drag forces effectively reduce the charged grain velocity. 
Up-to-date sputtering yields for graphite and silicate (olivine) grains have been derived 
using the code TRIM, for which we provide analytic fits. After training our method 
on a homogeneous density case, we analyze the grain motion and sputtering in the 
IGM density field as derived from a ACDM cosmological simulation at z = 3.27. We 
found that only large (a > 0.1-/im) grains can travel up to considerable distances (few 
X 100 kpc physical) before being stopped. Resulting metallicities show a well defined 
trend with overdensity S. The maximum metallicities are reached for 10 < (5 < 100 
(corresponding to systems, in QSO absorption spectra, with 14.5 < logA^(IIi) < 16). 
However the distribution of sputtered metals is very inhomogeneous, with only a small 
fraction of the IGM volume polluted by dust sputtering (filling factors of 18 per cent 
for Si and 6 per cent for C). For the adopted size distribution, grains are never com- 
pletely destroyed; nevertheless, the extinction and gas photo-electric heating effects 
due to this population of intergalactic grains are well below current detection limits. 
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1 INTRODUCTION 



In recent years, the possible presence of dust in the 
IGM has received new atte ntion. Before the observa- 
tions of CMB an isotropies ide Bernardis fc et alJ l200Ct 
ISpereel et a'i]l2003h strengthened the hypothesis of a flat, 
A-dominated universe, an explanation afternative to cos- 
mology was proposed for the dimming of distant (z ~ 0.5) 
Type la SNe jRiess et alJll998l : IPerlmutter et al.lll999^ : ob- 
scuration from intergaiactic dust may have been responsi- 
ble for the observations, provided that dust grains h ad a 
grey ex tinction law (as for elongated or large grains: Aeuirrd 
Il999bllah and a smoot h distribution tra cing the iower den- 
sity intergalactic gas llCroft et alJl2000ri . Although now the 
cosmofogic origin for the dimming is preferred, dust is stiff 
considered to be an important component of the extragafac- 
tic medium. Besides the most commonfy studied extinction 
effects, dust may affect, for instance, the thermaf bafance 
of the gas, being an effective heating agent in fow density 
environments (through photoefect ric heating; iNath et af J 
I1999I : llnoue fc KamavalEooal |2004|) and providing a source 



for coofing in the denser, higher temperature In tracfuster 
Medium (via IR emission; iMontier fc Giardlbooil . 

Furthermore, dust is studied for its own nature. About 
haif of the metafs in the Mifky Way and in other focaf gafax- 
ies is focked up in dust grains llWhittet|[T992: James et afJ 
l2002h . Though originating in stars, metafs are afso present 
in the fow density IGM, far from production sites. Metaf 
fines are now routinefy found associated to hydrogen in the 
z = 3 LyQ forest for cofumn densities as fow as N(Hl) 
< 10^"*'^, correspondi ng to gas onfy 10 times more dense 
than t he cosmic me an jCowie fc SongaifaliggllEffison et afJ 
I2OOOI: ISchavdIiool . The most popufar mechanism to ex- 
pfain the metaf poffution of the IGM is galactic winds 
from supernovae expfosions: metaf enriched gas is bfown out 
from (proto-)gafaxies into an IGM with pristine conditions 
jMadau et af 1120011 lAguirre et afj|200ldfbl) As an afterna- 
tive. iFereare^^t F lIlQQlf l and lAguirre et all (|2001a"") have 
proposed that metafs coufd be expeffed from gafaxics as 
radiation-pressure driven grains. Grains then refease metals 
in the IGM through sputtering. Not onfy this mechanism 
can in principfe provide for the same fevefs of metafficity as 
the galactic winds hypothesis, but the attractive feature of 
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this model is that, unlike winds, enrichment by dust would 
not impact the thermal/structural properties of the IGM, 
as no shock waves are involved in the transport. In addi- 
tion, metals transported via this process are associated with 
cool gas, alleviating the need to explain how this can be 
achieved if metals are transported by the hot gas produced 
by SN shocks. 

In spite of such interesting features several questions 
await a quantitative answer: how far can dust grains travel 
in the IGM? What is the dominant sputtering mechanism? 
How much dust is needed to produce the observed metallic- 
ity? These are the basic answers we are aiming at obtaining 
in the present work. 

A few authors have studied the ejection of dust grains 
from galaxies due to the radiation pressure from the galaxy's 
starlight. References to these works are given in § |21 Al- 
though no estimate is made about what fraction of the to- 
tal dust content is ejected, most works agree that escaping 
grains will be able to overcome the gravitational attraction 
reaching large velocities, v >100 km s~^. When colliding 
with gas particles, the kinetic energy of the impact is com- 
parable to the thermal energy of a gas with temperature 
T «10''-K. Thus, non-thermal sputtering may be a viable 
mechanism to erode grains and deposit their constituents 
(metals), even in the low temperature (and density) gas. In 
this paper, we will assume that a fraction of galactic dust 
grains will be able to go beyond the virial radius of a galaxy, 
with velocities and radii as inferred from literature. From 
this point onward the grain is slowed down by drag forces: 
we include both coUisional drag and Coulomb drag (due to 
the ionised nature of the medium and the charge a grain 
attains when exposed to the metagalactic UV background). 
During the grain motion, metals are deposited in the IGM 
as a result of sputtering (we consider both thermal and non- 
thermal processes). O ur approach is different from that of 
lAguirre et alj i2001al l^ , as they assume that dust grains can 
reach the equilibrium point between gravitation and radia- 
tion pressure and they include thermal sputtering only. 

The paper is organised as follows: in §|21we summarise 
the literature results on dust ejection from galaxies, which 
we use as initial conditions in our computation; §|3describe 
the physics we have adopted to describe the grain motion, 
charging and sputtering. Results are presented and discussed 
in § 01 for the ideal case of a homogeneous density field; 
the method is then applied to a more realistic simulated 
cosmological density field in § |S] In § |S| we will compare 
our results with previous works. Finally, we summarise our 
results in § . 



2 DUST EJECTION FROM GALAXIES 

Following the motion of dust grains from their formation 
sites in a galactic disk to the outer reaches of the halo is a 

^ Throughout this paper we will assume a fiat universe with to- 
tal matter, vacuum, and baryonic densities in units of the criti- 
cal density of £7m = 0.3, Oa = 0.7, and Q,t,h^ = 0.028, respec- 
tively, and a Hubble constant of = 100 ft km s~^ Mpc~^, with 
h = 0.7. The parameters defining the linear dark matter power 
spect rum are erg = 0.9, n = 0.93, dn/dlnk = CSpergcl ct aU 
|2003^ 



complex task. It involves the evaluation of: radiation pres- 
sure a nd other forces resulting from anisotropic radiation 
fields llWeingartner Sz Drainel[2001lJl . possibly taking into 
account the opacity of the dusty disk IIDavies et al. 1993); 
disk and halo gravity; gas drag (and the grain charge; see 
next Section); sputtering rates; magnetic forces; geometry 
and physical condition in a multiphase disk and halo gas. 

A few authors have studied the dust expulsion, includ- 
ing most (but not all) of the relevant processes. We sum- 
marise here their findings. 

(i) Depending on their size and composition, some dust 
grains can escape the halo of the galaxy in a rela- 
tively short time (a few hundreds Myr). These escaping 
grai ns attain termin a l velocit i es of about 1 00-1000 k m s~^ 
( Fe rrara et al.l Il99ll: iFerraral Il997l : IShustov fc Vibe) Il995l : 
I Simonser^^Hmm^stadl^ggl) ) . 

(ii) Escaping grains are relatively large, with sizes in the 
range 0.05-0.2 /im. Larger grains are too heavy to escape the 
gravitational well, while smaller grains offer smaller efficien- 
cies to radiatio n pressure and do not travel far from thei r 
formation site dShustov fc Vibdll995l: iDavies et al.lll998f) . 
Furthermore, smaller and slower grains could be more effec- 
tively eroded by sputte ring, as they spend more time in the 
hot halo environment jFerrara et al.lll99lt IShustov fc Vibl 
Il995h . Instead, velocitie s for escaping grains do not d epend 
much on the grain size llSimonsen fc Hannesta 

(iii) Being heavier and having smaller radiation pres- 
sure efficiencies, silicate grains reach smaller velocities than 
graphite grains and may be underrepresented among the 
escaping grains iBarsella et al.l Il989t iFerrara et alJ ll99Cl . 
1991). 

Unfortunately, none of these works provide a statistic 
for the properties and amount of escaping grains. Apart from 
a few illustrative cases, here we will simply assume that sizes 
are in the range 0.05< a [fim] <0.2 and velocities in the 
range 100< v [km s~^] <1000, with all values of a and v 
equally represented. Typically, dust grains in a galaxy are 
believed to obey a power law size distribution, n(a) oc a~^'^, 
and to have a ratio 1:1 for the relative proportion of graphite 
and silicate grains (i.e. the MRN model for Milky Wa y dust; 
iMathis. Rumpl fc Nordsiekll977HDraine fc Leeil984) . For a 
given dust mass, the MRN model sets the number of grains 
for each radius bin, with smaller grains being obviously more 
represented. By adopting a flat distribution for the sizes of 
the grains that eventually escape in the IGM (i.e. with the 
same number of large and small grains), we implicitly as- 
sume that a good fraction of the dust mass is trapped inside 
the galaxy. Even allowing all large grains to be ejected, at 
maximum only about 10 per cent of the dust mass can travel 
to the IGM. Thus, we use 10 per cent as an upper limit to 
the fraction of a galaxy's dust mass that can be ejected. The 
limit has been derived for the MRN graphite and silicate 
abundance ratio 1:1, that we adopt also for ejected grains. 
As literature results suggest a lower efficiency for ejection of 
silicate grains, this ratio may constitute an upper limit to 
the presence of silicate grains in the IGM, if the assumption 
of MRN-like dust inside a galaxy is valid. 

Unfortunately, size and material distributions for dust 
in external galaxies are quite uncertain, especially at high 
redshift. In the early universe, du st can be efficiently pro - 
duced only by Type II supernovae llTodini fc Ferrarall20oih . 
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For pristine conditions of the parental cloud, dust is pro- 
duced in pair-instability supernovae from massive (140-260 
Mq) stellar progenitors: models predict in this case a pref- 
erential formation of silicate grains (iNozawa et alj l2003t 
ISchneider et alj|2004) although considerable masses of car- 
bon can be locked up in dust i f the CO molecule fo rma- 
tion is duly taken into account dSchneider et aljl2004t) . As 
the gas metallicity ri ses above a critical value (> 10~^ Zq; 
ISchneider et alJl2003D the formation of very massive star is 
suppressed and dust forms in the ejecta of core-collapse su- 
pernovae: in this case, similar masses of carbon and silicate 
grains can form Ijodini & Forrara 20^. For both kind of 
objects, the cumulative size distribution over all the materi- 
als formed in the ejecta can be roughly described by a power 
law; in pa rticular, large grains tend to follow the MRN dis- 
tribution jNozawa et al.ll2003h . Thus, we believe the upper 
limits derived above for the amount of dust ejected into the 
IGM are reasonable estimates. 



3 PHYSICS OF DUST GRAINS IN THE IGM 

In this Section we present the basic physical ingredients we 
have adopted to describe the motion, charging and sput- 
tering processes for dust in the IGM. We will consider a 
gas with pristine composition (76 per cent H, 24 per cent 
He in mass). Before considering a cosmological density field, 
we will show some illustrative cases for the mean density 
at z = 3. For the adopted cosmology, the mean baryonic 
density at z = 3 is pb = 2.3 x 10~^^ g cm~^ (equivalent 
to riH = 1.0 X 10^^ cm~^). As for the gas temperature, we 
have adopted T = 2 x 10''-K, a typical value d erived for the 
low d ensity IGM from QSOs absorption lines iSchave et alJ 
|200[|) . As we allow for different ionization fractions, six dif- 
ferent kind of particles may be present in the IGM: neutral 
and ionised hydrogen (Hi, Hll), neutral, single and double 
ionised helium (He I, Hell, He ill) and electrons. 

We will find useful to express the grain velocity v 
through the atomic speed ratio 



ruiv 
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where rrii is the mass of each particle constituting the gas. 
It is SH=SHc/2=-y/mH/m.oSo ~42.8sc. We saw in §|3that 
escaping grains achieve a certain terminal velocity as a bal- 
ance between radiation pressure and the other forces op- 
posing the motion. That velocity will be considered as the 
initial velocity of the grain in the IGM, and we will follow 
the motion assuming that dust is slowed down because of gas 
drag ('i) l3.1ll . For simplicity, the gas drag is the only force in- 
cluded in our study. This is clearly an approximation when 
studying the motion of a grain in a inhomogeneous density 
field (§ EJ. In that case we should also consider the path 
deviations due to gravitational attraction by denser regions; 
if the encountered overdensities are forming stars, radiation 
pressure an d other forces resulting from anisotropic radi- 
ation fields llWeingartner fc Draing l2001b^ should be also 
included. However, the approximation is not too crude: at 
the virial radius of a z = 3 galaxy of median mass (§01 
and §|5J|, the gravitational and radiation pressure force (as- 
suming a galactic mass-to-light ratio) on a a—O.l-fim grain 
nearly compensate each other, being of order 10~^^ dyne. 



The drag force due to the IGM gas is of the same order (see 
Fig.0 and dominates as the grain moves from its initial po- 
sition, since the drag force on high velocity grains decreases 
less than the distance squared. 

Grains in the IGM are charged, because of collisions 
with ions and because of the photoionization due to the 
metagalactic UV background. The charge influences both 
the drag and the sputtering efficiency. The charging pro- 
cesses we have considered and the typical charges of a grain 
in the IGM are described in § 13.21 We will denote with U 
the potential on the grain surface {U = Ze/a, where Z is 
the grain charge in units of the electric elementary positive 
charge e; Z — 694.5f/[Volt]a[^m]). It is helpful to define the 
reduced potential 

^ KT 

Charged grains could be deviated from their path by the 
Lorentz force, if moving through a magnetic field. For sim- 
plicity we neglect this force here. If the primordial IGM 
magnetic fields are l ow as predicted in c osmological simula- 
tions {B « 10"^^ G: lGnedin et alJl200nh . the Larmor radius 
would be of order 10 Gpc for the grain properties consid- 
ered in this work (S I3.2.3|l . Being this much larger than the 
typical length travelled by a grain in a Hubble time, the 
Lorentz force could indeed be neglected. Unfortunately, it 
is not clear if efficient large scale amplification mechanisms 
are present. Some observations seem to indicate intergalac- 
tic B-fields as high as 0.1 nC (iVall oe 2004, and references 
therein), in which case charged grains would preferentially 
move along field lines. However, the inclusion of the Lorentz 
force is still prevented by the large uncertainties on the field 
direction, strength and structure. 

The sputtering process is described in § 13.31 Because 
of collisions with gas particles, atoms may be knocked off a 
dust grain and released in the IGM. In a low temperature 
medium, the most relevant mechanism is non-thermal sput- 
tering, in which the energy transfered to the grain's atoms is 
the kinetic energy of gas particles impacting at a supersonic 
grain speed v. Since we consider grains made of graphite 
and silicate (olivine), sputtering could be able to pollute the 
IGM with carbon, iron, magnesium, silicon and oxygen. 

Finally, § 13.41 shortly describes how the physical pre- 
scription are implemented in our numerical code. 



3.1 Drag Forces 

A charged dust grain moving through a gas experiences a 
drag force due to the direct collision with the particles in 
the gas and to Coulomb interactions with the ions. Under 
the hypotheses that the collisions between grains and ions 
are elastic and that the grains are much smaller than the 
mean free path of the p articles in the gas, we can follow 
IPraine fc Saloeteil (|l97flh and write 



- drag 



2Tva^KT 



1^ rii [Go(sO + ln(lA/^,l)G2(sO] | , (1) 

where the summation is over the six kind of particles we 
have considered for the gas, Zi is the charge of each particle 
(in units of e; Zi = — 1 for electrons) and 
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Figure 1. Drag force on a spherical dust grain of radius a = 0.1- 
/im moving at velocity v through a gas with T = 2 X 10*-K and 
2 = 3 mean cosmic density. The gas is composed of Hll, Hell 
and electrons. The dotted lines is the force component due to 
collisional drag. The dashed line is the component due to Coulomb 
drag for a grain potential (7=20- Volts. The solid lines represent 
the total drag. 
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The f irst part of Eqn. l[T| is the term due to the collisional 
drag iBaines. Williams fc Asebiomdll965ll . ignoring a neg- 
ligible contribution due to Coulo rab focusing. The second 
part describes the Coulomb drag iSpitzeilll962h . For Go(s) 
and G2{s), whose exact formulae depend on the error func- 
tio n, we use the simple ana lytical approximations derived 
bv lDraine fc Saloeteil Jl979ft . 



Go(s) ; 



30F 



1/2 



which provide an accuracy within 1 and 10 per cent, respec- 
tively, for < s < CX3. 

In Fig.Qwe show the drag force acting on a dust grain 
of radius a=0.1-/xm as it moves with velocity v through a 
gas of mean z — 3 cosmic density and T = 2 x W*-K. 
We have assumed that the IGM is composed of Hll, Hen 
and electrons. If the grain is not charged (or the gas not 
ionised) only the collisional drag is present (dotted line). 
For a supersonic dust grain (sj ^1) the collisional term in 
Eq.El reduces to llShulllll97^ 



-fcoll 



(2) 



As it can be seen in Fig. Q the supersonic regime is reached 
at low velocities: for H and He we have v 3> 18.1 km s~^ 
and V ^ 9.0 km s~^, respectively (for the gas temperature 
adopted), and -Fcoii oc beyond that velocity. At high veloc- 
ities, collisions with helium ions constitute about 30 per cent 
of the collisional drag. Free electrons, instead, do not carry 
enough momentum to contribute more than a few percent 
to the drag force, at any velocity. Thus the collisional drag 
is nearly independent of the ionization state of the IGM. 



For a charged grain moving through a plasma, one needs 
to consider also Coulomb drag. This is shown in Fig. 
(dashed line) for a grain with a typical potential U = 20- 
Volts (see S 13.2.31 the Coulomb drag does not depend on the 
sign of the grain charge). For small velocities, the Coulomb 
drag increases linearly with v. For supersonic velocities, in- 
stead, it decreases as . The first maximum at v ~ 15 km 
s~^ occurs when the dust grain has a velocity similar to the 
thermal velocity of H and He (sh ~ 1 and sho ~ 1). For 
larger w, the contribution of H and He to the Coulomb drag 
decreases until the electron component become dominant. A 
second local maximum is reached when So = 1 (w — 680 km 
s~^), after which the drag decreases again. 

It is interesting to note that the Coulomb drag de- 
creases with increasing gas temperature. This is because 
of the increasing number of gas particles moving faster 
than the grain, which do not contribute to the net drag 
iNorthrop fc Birmingha ml ll990ri . The collisional drag, in- 
stead, depends on y/T for low velocity and is independent 
of T in the supersonic case. 



3.2 Grain Charge 

Dust grains immersed in a hot gas and subject to a UV 
background attain an electrical charge because of: i) colli- 
sions with electrons (which tend to make the charge more 
negative); ii) collisions with positive ions, and iii) photoejec- 
tion of electrons from the grain by absorption of UV photons 
(which tend to make the charge more positive) . If we denote 
with Jci the charging rate due to collisions with ion i (i.e. 
the number of charges, in units of e, captured by a grain 
colliding with particle i per unit time) and with Jpc the rate 
for photoelectric charging, the grain charge at equilibrium 
can be found by imposing 



0. 



We ignore here charge quantization JPraine fc SutiJl98itl . 
This is a safe assumption for the large grains adopted in this 
work and for the large potentials they attain (see § 13.2.31 . 



3.2.1 Collisional Charging 

We foUowed lSMiil lll978f) to derive the charging rates on a 
a spherical dust grain that moves relative to the gas. The 
charging rate due to collisions with particle i is^ 

KT 1 



Jci = TTa^ni{zi£_i + Si) (- \ 

\ Inmi J Si 

+ s? - [erf (si - So) + erf (si + So)\ (3) 
+ (Si + So) exp[-(si - Sof \ -f (si - So) exp[-(s, -I- So)^]| 

with erf the error function, the sticking coefficient of the 
particle i, Si the secondary electron emission coefficient and 



^ Electrostatic polarization of a du st grain by the electr ic field 
of an approaching charged particle iPraine fc Sutinlll987l) is not 
included in Eqn. |3. Its effects are small for the typical charges 
described in S 13.2.31 
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1 ^Ziip Zi(t) > 

the minimal atomic speed ratio for the approaching par- 
ticle to win electrostatic repulsion. In the limit Si 0, 
Eqn. lEl redu ces to the classical result for a static grain 
JSpitzeJligfl) . If the grain moves relative to the gas, the 
frequency of collisions with ions and electrons is altered. For 
Si 3> 1, Eqn. Q tends to 



1/2- z,<j>. 



(4) 



As for the sticking coefficients, we assume that half 
of the colliding electrons are captured by the grain (^e = 
0.5), the other half being scattered; and that all the ions 
will neutralise as they arrive on the surface, thus sharing 
their positive charge with th e grain (Ci = 1: iDraindHoTa 
IWeingartner &: Drainel *2001c'). For the coefficients for the 
secondary emission of electrons, Sj, we used the empirical 
expressions given bv lDraine fc Salpeteil (|l979). The coeffi- 
cients depend on the energy of the impact, {Eq) + 
wher e (Eo) is the mean therma l energy of the impinging par- 
ticle dPraine fc Salpeteilll979|) . Because of the small sput- 
tering rates fij 13. 3t . we neglec t the charge taken away from 
the grain by sputtered atoms jPraine fc SalpetejEgTSl) . 



3.2.2 Photoelectric Charging 

The photoelectric charging rate can be written as 



'^pct 



?abs(/ii', a)Y{hv, a) ^^'^'^ du, 
hv 



(5) 



where Jv is the mean specific intensity of the UV back- 
ground, Qahsihv, a) is the absorption efficiency of a dust 
grain of radius a and Y{hv,a) is the photoelectric yield, i.e. 
the probability that an electron is ejected when a photon 
absorption occurs. The lower limit of integration is the pho- 
toelectric threshold frequency, I'pct, for which we use 



W + eU [/ > 0, 



W 



U <Q, 



W being the workfunction, i.e. the ionization potential of 
a neutral bulk material. Thus, if the grain has negative 
charge, an electron is photoejected as soon as a photon 
with energy hu > W is absorbed. If the grain has a pos- 
itive charge, the photoejected electron needs to have a ki- 
netic energy large enough to escape the attraction of the 
charged g:rain, therefore ionizing photon s must have hv > 
W + et/ . IWeingartner fc Draiiij feOOldl derive the photo- 
electric threshold for spherical grains taking into account 
geometric effects and polarization. They obtain values for 
Upet that differ from those used here by additional terms of 
order e^/a and dependent on . For the large grains used 
in this work, those terms are negligible compared to W and 
we omit them here. We use W — 8- eV and W = 4.4-eV for 
graph ite and silicates, respectively iWeingartner fc Draind 
12001 d) . 

Following IWeingartner fc Drainel ll200ld) . we write the 
photoelectric yield as 



-23 
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Figure 2. The z=3 UV background of lBianchi et alj <200ll) in- 
cluding the contribution of galaxies and QSOs (solid line). The 
vertical bars show the workfunction adopted for silicates and 
graphite. The dotted line is the QSOs contribution to the back- 
ground. 



with yo the photoelectric yield of the bulk material and O 
hi' — W. For graphite and silicates we use, respectively 



o.oo9(e/w)^ 

1 -H0.037(e/Vl/)5 



o.5(e/iv) 
i + 5{e/w)' 



In small particles, the photoelectric yield is enhanced with 
respect to bulk materials. For spherical grains of radius a, 
the enhancement factor is well approximated by 



yi{hv,a) = 



/32 - 2/3 2 - 26-/3 



mm. 

l200ldl . 



with j3 — a/ la and a = a/la + a/l^, la and 1^ being t he pho- 
ton attenuation length and the electron escape length Draine 
We adopted le = lOA jWeingartner fc Draine 
Using tabulated values for the optical proper- 
ties of graphite and 'smoot hed astronomical silicates'^ 
llWeingartner fc Draindf2001a^ . we derived the photon at- 
tenuation length with la = A/[47rIm(m)], where A is the ra- 
diation wavelength and m(A) the complex refractive index. 
A proper weighted mean h as been used to take into accou nt 
the anisotropy of graphite (IWeingartner fc Draindl200lJ) . 

Laboratory measures on bulk materials have shown that 
the distribution of kinetic energy for photoejected electrons 
drop to zero at E = and E = hv — W and peaks at 
interme diate energies iDraindll978ll . lWeingartner fc Draind 
J200ldl adopted a parabolic energy distribution and derived 
the probability for electron escape to infinity, 2/2. Again ig- 
noring all terms of order e^/a and dependent on a~^, it is 



y2{hv,a, Z) 



■ ) hv-W > 



C/>0, 
(7 < 0. 



Finally, we have used tabulated values of Q^]^s(hv,a) 
for spherical grains of graphite and 'smoot hed astronom- 
ical silicates' 'Weingar tner fc Draind lioOlah . We did not 
consider photodetachment of electrons in the energy lev- 
els above the valence band of negatively charged grains 
l|Weingartner fc Drainel2001a^ , which can be ignored for the 
grain charges and radii considered here. 



Y(hv,Z,a) = y2(hv,a,Z)imn[yo{Q)yi{hv,a), 1], 



Available at |http : //www ■ astro ■ prlnceton ■ edu/~draine | 
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Figure 3. Equilibrium potentials for graphite and silicate dust grains exposed to the 2 = 3 UVB, as a function of the gas overdensity. 
In the left panel we show the equilibrium potentials for static (ti = 0) grains of radii 0.01, 0.05, 0.1 and 0.2— /im (in both panels, the 
upper bundle of curves refers to graphite grains, the lower to silicates). We have assumed a gas made of Hll, Hell and electrons, with 
temperature T = 2 X 10*-K. In the right panel we show the effects of variation in grain velocity and gas temperature for a grain with 
radius a=0.1-/im. For each grain material, the solid line refer to the same conditions as in the left panel, the dashed line to a grain 
moving with v = 500 km s~^ in a T = 2 X 10*-K gas, the dotted line to a static grain in a T = 2 X 10^-K gas. 



3.2.3 Typical grain charges in the IGM 

We adopt in this work the ultra violet background (UVB) of 
iBianchi. Cristiani fc KTiiil J200 J) . The UVB, shown in Fig.lH 
includes contributions from both QSOs and galaxies, assum- 
ing that a fraction /esc ~ 0.1 of the Hi- ionising UV Ught 
can escape internal absorption. For any observer's redshift, 
the UVB shows two breaks at frequencies v > 13.6-eV and 

V > 54.4-eV in the observer's rest frame, due to the absorp- 
tions of Hi and Hen Lyman continuum photons, respec- 
tively, from the residual neutral gas in the IGM. While the 
QSO contribution (dotted line) can be fitted between the 
breaks by a power law (mainly because of the assumption 
on the intrinsic QSO's spectrum), the Hi- ionizing UVB in- 
cluding galaxies is only roughly described by cx u~^'^ for 

V < 54.4-eV. The contribution to the photoelectric charg- 
ing rates of frequencies v > 54.4-eV is negligible, therefore 
we use fmax ~ 54.4-eV in Eqn. @. The UVB is computed 
integrating light coming from all objects up to a maximum 
redshift Zmax- In iBianchi et al. ( 2001) we adopted Zmax ~ 5. 
This assumption has no effect on H l-ionising photons: be- 
cause of IGM absorption only local sources contribute to the 
calculation. For radiation at lower frequencies not able to 
ionise Hi, the UVB depends on the redshift assumed for the 
first sources, mainly because of the near ly constant galaxy 
emissivity at 2: > 2 jBianchi et al.ll200ll) . However, as long 
as we use the UVB at 2; < 4, the photoemission rate in 
Eqn. ||^ is only slightly affected by our choice for Zmax- 

In Fig.|3we show the equilibrium potentials for graphite 
and silicate grains exposed to the z — 3 UVB, as a function 
of the gas overdensity. Again, the gas is composed of Hll, 
Hell and electrons and has T — 2 x 10*-K. In the left panel 



equilibrium potentials are shown for static {v = 0) grains of 
radii 0.01, 0.05, 0.1 and 0.2-/im (spanning the range of radii 
relevant to this work). Because of the different photoelectric 
yields and absorption efficiencies, graphite grains generally 
attain a larger (more positive) potential than silicate grains. 

In the low density IGM the grain potential reach a rel- 
atively high positive value (a few ten s of Vo lts), in agree- 
ment with what found bv iNath et"ai] (l99^. The equilib- 
rium potential is reached when the coUisional charging rate 
for electrons Jcc is balanced by the photoelectric charging 
rate Jpc, the contribution of positively charged ions being 
greatly reduced by Coulomb repulsion {Jd ^ Jpc). While 
the coUisional charging rate depends on the particle density, 
the photoelectric charging does not. Thus, for low values of 
We, a high Coulomb attraction is needed to increase the col- 
lisional cross section, in order to have Jce ~ Jpe- Charges for 
the grains studied here are always much lower than the max- 
imum limit over which fi eld emission of positive ions occurs 
JPraine fc Salpetei|[l979l) . 

At higher density, the contribution of collisions to charg- 
ing becomes progressively more important and the charge 
decreases. This partly explains why the slope of the equi- 
librium potential of a grain of given radius changes at 
higher densities (see Fig. |^. A second reason for this is 
that i/pet goes across the Lyman discontinuity in the UVB. 
For very high densities, only coUisional charging rates are 
important. Neglecting secondary emission, the classical re- 
sult is recovered: grains have negative potent ial because of 
the higher collisi on rate of electron in a gas (ISpitzedllQTi : 
IOsterbrocldll989ll . For the plasma composition and sticking 
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coefficients adopted here, the equifibrium would be reached 
for (j) ^ -2.0. 

For a grain that moves at a supersonic speed, the col- 
lisional charging rate tends to Eqn. ||1J and becomes larger 
than that in the static case. When Si 3> 1 but still Se ^ 1, 
Jci gives a constant contribution to the positive charging 
rate, with a small dependence on the grain potential. When 
also Se ^ 1, Jce and Jd almost compensate, apart from 
the contribution of gas particles impacting from directions 
normal to the grain motion, which results in a net negative 
charging rate proportional to (p/Se- Because of this, mov- 
ing grains always have higher (more positive) equilibrium 
potentials. In Fig. |3| (right panel) we show the equilibrium 
potentials for a a=0.1-/im grain moving in a T = 2 x 10*-K 
gas Sit V — 500 km s~^. For the velocities explored in this 
work, differences with the static case are never larger than 
about 5- Volts. 

In Fig. 121 (right panel) we also show the change in the 
equilibrium potentials of a a=0.1-/im static grain when the 
temperature raises to T = 2 x 10^-K gas. At higher tem- 
peratures gas particles move faster and the Coulomb cross 
section reduces. At low densities, a reduced Jca has to com- 
pete with a Jpc unaltered by temperature changes and the 
equilibrium potential increases. According to the classical 
result, at high densities the charge would reach a high nega- 
tive value. This is not the case when the secondary emission 
of electrons is included. The secondary emission of electrons 
decreases the negative charging rate due to electrons Jcc 
and increases the positive charging rates for ions Jd. As 
a results, equilibrium potentials are higher when this pro- 
cess is included. The effect can be larger for higher tem- 
peratures (or for supersonic grains). In Fig. |3 the positive 
equilibrium potential attained in the high density gas for 
T = 2 X 10^-K is due to secondary emission. If Sc,5i—0, a 
static grain would have an equilibrium potential U — —34.5- 
Volts, while it is 4.2 and 9. 8- Volts for graphite and silicates, 
respectively, when the secondary emission is accounted for 
(silicates have higher secondary emission for impinging elec- 
trons). Thus, secondary emission prevents dust grains in a 
hot overdense gas to reach high negative charges. The lim- 
iting negativ e potential below which field emission of elec- 
trons occurs dWeingartner Sz Drainel200ld) is never reached 
in this work. In the T = 2 x lO^-K gas, secondary emission 
raises the equilibrium potential by less then a few Volts, and 
its effect can be mimicked by reducing the electron sticking 
coefficient to ~ 0.4 (in most cases, the secondary emission 
due to colliding electrons is more important than that due 
to positively charged ions). 

As the parameters describing the galactic contribution 
to the UVB (namely, the star formation history and the /esc 
fraction) are quite uncertain, especially at high-z, we also 
computed charges when only the QSO contribution (dotted 
line in Fig.|2J is taken into account. Because of the reduced 
UV fiux, the equilibrium charges are smaller in this case. 
However, the diflerence is small (a couple of Volts) and it 
does not modify significantly the results we present in the 
rest of the paper. 

3.3 Grain Sputtering 

The efficiency of the sputtering process is given by the sput- 
tering yield, Yi{E, 9), i.e. the number of atoms or molecules 



of the target material that are sputtered in each colli- 
sion with a projectile atom of material i. For each tar- 
get/projectile combination, the sputtering yield depends on 
the impact energy E and on the angle 6 (relative to the 
surface normal; 6 = Q for normal impact) of the impacting 
particle. Our targets are dust grains made of graphite or sil- 
icates, while we will consider neutral and ionised hydrogen 
and helium atoms as projectiles. The sputtering yield does 
not depend explicitly on the charge of the projectile (ions 



rapidly neutralise when approaching a solid surface; 
although the energy of the impact E depends on Coulomb 
repulsion / attraction. 

From Yi we can derive the sputtering rate, i.e. the num- 
ber of atoms that are sputtered off each grain per unit time, 
dN/dt. For a charged (non-rotating) spherical dust grain of 
radius a that moves with velocity v through a Maxwellian 
gas, we can write (Draino & Salpoter 1979) 



dN 
~dt 



g.{e.,4>){Y,[E = {e-z.<i,)KT])ode, (6) 



where again the summation is over all relevant projectile 
particles of species i, emin = max[0, and {Yi{E))e is the 
angle averaged sputtering yield, given by 



r/2 

{Yi{E))e = 2 j Yi{E, 6) sine cos Ode. 



(7) 



The function gi is given by 



5i(e,(/>) = 



exp - ej rr f Zi(t> 



^ (l - ^) sinh(2x/^s0- 



For supersonic grains (si oo) the function Qi tends to 



5.(e,0) = (l-^^)5(e 



and the rate for non-thermal sputtering is obtained, 

dN 
~dt ~ 



\ynAl~'^^){Y,[\m,v^~z,eU])e, (8) 



where the term within round brackets is a co rrection to th e 
grain cross section due to Coulomb focusing JSpitzeilll97'3l . 

A few analytical models for the sputtering yield of 
astro nomical dust mate r ials can be found i n the litera- 
ture JPraine fc SalpeteJ [l979l: iTielens et ahl Il994^ . As a 
comparison, we refer h ere to the sputtering yields of 
iDraine fc Salpeteil il979|) for graphi te and for olivine (a typi- 
cal silicate, FeMgSiOA: lDrainei l995'). Normal incidence sput- 
tering yields are shown in Fig. 0] together wit h a few exper- 
imental results which are taken from lDraiii3 lll995l) . Above 
a certain threshold energy, the yield increases steeply, reach 
a maximum and then decreases slowly for higher impact en- 
ergies. The normalization of Yi{E,9 = 0°) was chosen to fit 
the available data. Because of the lack of experimental data 
for low impact energies, which is the most important in as- 
trophysical applications, the sputtering yield in proximity of 
the threshold is poorly constrained. The dependence on the 
collision angle is uncertain as well, and it is usually assumed 
that Yi(E,9) ~ l/cosS. Under this assumption, it follows 
from Eqn. ||7J that the angle averaged sputtering yield is 
twic e the norma l one. 

iField et al.l I^1997^ and iMav et alJ (|200nl: see also 
iFlower et al.l 19961 : lJurac et alJ Il998l) derived sputtering 
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Figure 4. Normal incidence (8 = 0°) sputtering yields for H and He on 
energy. Solid curves are theoretical estimates from lPraine fc SaloeteJ il9' 



iraphite and sili cates (ol i vine) , as a function of the impact 
graphite) and iDrain j ilQQfil . olivine). Small dots are our 



resul ts using the T RIM 2003 code, together with the fits of Eq. J5J (dotted lines; see text for details). Squares are experimental measures 
from lPrainj ll99fi) : closed symbols are for He, open symbols for H. As there arc no measures for sputtering on silicates, datapoints in 
the right panel are for a similar material. Si O2 iDrainclQQM . 



yields using the M onte Carlo code for the TRansport of Ions 
in Matter (TRIM; IZieder et aT]ll985^ . Briefly, TRIM simu- 
lates the bombardment of a plane-parallel target by projec- 
tiles of given kinetic energy and impact angle with respect 
to the surface. The target is simulated as amorphous (i.e. 
given the density p, the position of each atom inside the 
target is random). The depth inside the material at which 
the projectile hits a target atom is derived from interatomic 
potentials dependent on the colliding elements. If an atom 
in the target is given an energy larger than the Displace- 
ment Energy En, it will be able to escape from its position 
(the lattice site) and it will lose to the lattice an amount 
of energy given by Eb, the Bulk Binding Energy. Recoil- 
ing atoms may collide themselves with other target atoms. 
TRIM follows the whole collision cascade. Finally, an atom 
in the proximity of the target surface is counted as sput- 
tered if the component of its kinetic energy normal to the 
surface is larger than the Surface Binding Energy, Es- Sev- 
eral projectiles are needed to derive a mean sputtering yield. 
In order to have an estim ate alternative to the widely used 
iDraine fc SaloeteJ ^1979^ yields, we have used the 2003 ver- 
sion of TRIM* to derive the sputtering yield for graphite 
(p= 2.266 g cm"^, E p = 25 eV, Eg = 3.0 eV, Es = 7.41 
eV; iField et al.lll997D and olivine (p=3;843_scni;^ Ed = 
50 eV, Eb = 9.7 eV, Es = 5.64 eV; lMav et aPbOOOD . 

We ran TRIM to have Yi{E,e) for several values of E 
and d. Our results for sputte ring of olivine by He are very 
similar to those published bv lMav et al.l (l200(i) . The sput- 
tering yields we obtain for He into graphite are at least a fac- 



* The TRIM program is part of the SRIM package, which can be 
downloaded at |http : / /www . sr im ■ org| 



tor of two larger than those reported bv lField et ai] Jl997D . 
probably because they used an earlier version of the code. 
None of those author gives sputtering yields with H as a pro- 
jectile. In Fig. 13] we show the sputtering yields for 9 = 0°. 
Above a certain threshold energy, yields increase steeply, 
reach a maximum and decrease at a slower rate for higher 
energies (because projectiles implant deeper into the target). 
At high energies, the behavior o f the TRIM results is sim - 
ilar to the analytical models of iDraine fc SalpeteJ il979l) . 
although graphite yields for both H and He are higher than 
the laboratory results. In the proximity of the threshold, dif- 
ferences are larger. In general, the threshold energies of the 
analytical yields are lower than those derived with TRIM 
(with the exception of He into graphite). Such discrepan- 
cies may be the result of uncertainties in the estimate of 
the energy parameters {En, Eb an d Es) adop t ed in the 
TRIM calculations, as discussed by iField et al.l lll997l) . In 
particular, we found that a variation of 20 per cent in Es 
causes a variation of about 30 per cent in the sputtering yield 
at high energies, while the region of the threshold depends 
more on the choice of En'- by varying En of 20 per cent the 
threshold energy changes by 25 per cent We stress again the 
uncertainty in the derivation of the yield in the proximity of 
the threshold; the cross sections used by TRIM in comput- 
ing the collisions at low energies are based on extrapolations 
from experimental data which are only available at higher 
energies. 

The behavior of the sputtering yield for 6' > 0° de- 
pends on the energy E. For high E values, the yield in- 
creases with 8 faster than the usually assumed 1 / cos 8 de- 
pendence, while at energies lower than the energy corre- 
sponding to the maximum yield, Yi{E,6) is flatter and de- 
creases when the projectile approaches grazing incidence. 
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Table 1. Best-fit parameters to be used in Eq. |5]to calculate the angle-averaged sputtering yields of H and He on graphite and olivine. 
For olivine, parameters are given for each atomic species in the compound and for the total number of sputtered atoms. The parameters 
in the left panel refer to {Yi{E))0 as defined in Eq. J7J (the streaming yield adopted in this work). For completen ess, in the r ight panel we 
also give the parameters derived using Eq. IIUI to define the mean (the isotropic yield, preferred by^eld ct al. 1997] and^av et al. 2 00ol) . 
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As a result, {Yi{E))e > 2Yi{E,9 = 0°) at high energies and 
(Yi{E))e « YiiE ,9 = 0°) close to the threshold (see also 

alJ ll998h . TRIM calculations are 



the discussion in I.Turac et 
made using a plane-parallel infinite target. Sputtering yields 
may increase for spherical grains of th e dimension of th e 
mean penetration depth of a projectile jjurac et alJll99^ . 
As we mainly deal with large grains, we do not consider this 
radius-dependent enhancement. 

The data points computed with TRIM were fitted with 
the following function: 



Y{E) = fcexp 



E — Eth 



7 (in 



E 



En 



(9) 



Whil e the first term in Eqn. @ (introduced bv lMav et alJ 
I2OOOI) well describes the rapid increase of the sputtering yield 
in the proximity of the threshold energy Sth, the second 
term, centered on the energy value -Bmax, is needed to re- 
produce the maximum and the slow decrease of Y{E) at 
higher energies. The function provides a remarkably good 
description of both Yi{E, 9) (see, for example the TRIM data 
points for Y{E, 9 = 0°) and their fit in Fig.|3|l and the angle 
averaged yield {Yi{E))g of Eqn. (O, with most of the data 
points (~ 95 per cent) being within less than 10 per cent 
of the fitted function. The best-fit values of Ef^, £max and 
the parameters k, {3 and 7 can be found in Table (right 
panel) ^. For olivine, sputtering yields are slightly different 
for each of the elements in the compound. For simplicity, we 



^ The angl e- averaged yield o f Eq. J7J is referred to in lField et alj 
^1997^ and lMav et all i200(]t) as the streaming yield, appropriate 
for projectiles which streams in one dimension and collides with a 
non-rotating spherical grain. However, for a isotropic distribution 



use here the sputtering yield for the total number of atoms. 
Parameters for the fit of the total sputtering yield of olivine 
are also given in Tabled 

To show the differences between the analytical and the 
TRIM sputtering yields, we plot in Fig. |3 (right panel) the 
sputtering rates computed from Eqn. Q for a graphite or 
olivine grain moving through a gas composed of Hll, Hell 
and electrons. The grain has radius a — 0.1-fim and equilib- 
rium potential U = 20-eV. The gas has mean z = 3 cosmic 
density and temperature T = 2 x 10''-K. For both analyti- 
cal and TRIM sputtering yields, the low velocity sputtering 
rate is obviously dominated by collisions with the more mas- 
sive He atoms. The contribution of the more abundant H 
atoms becomes important for v > 100 — 200 km s~^ (where 
both He and H sputtering yields depend weakly on the im- 
pact energy). As said for the normal incidence case, the 
largest differences between TRIM and analytical sputter- 
ing rates are for graphite grains. However, the difference is 
small and the results shown in this paper do not change sig- 
nificantly when one of the yields is preferred over the other 
(see §BJ- For the mean conditions considered in this work. 



of projectiles it is more appropriate to use the isotropic yield, 
given by 



{y^iE))e 



r/2 



Y,{E,e)smede. 



(10) 



While lField et al.l <1997l) and lMav et~ai] <2000t) prefer the latter, 
we use in this work the former. As they pointed out, differences 
between the two means are not large. For completeness, we also 
give the best fit parameters for the isotropic yield in the right 
panel of Table [Tl 
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Figure 5. Sputtering rates for a graphite and olivine grain of radius a = 0.1-^m, moving with velocity v through a gas with mean z = 3 
cosmic density. The gas is assumed to be composed of Hll, Hen and electrons. In the left panel the sputtering rate s are shown for the 
TRIM results (graphite- solid line, oli vine- dotted-dashed line) and for the analytical model (graphite- dashed line, iDraine fc Salpeteil 
llQygt olivine- dotted line. IP rainellQQ^) . for a mean gas temperature T = 2 X Iff'-K and a grain equilibrium potential (7=20- Volts. In the 
right panel we show the TRIM 2003 sputtering rate for graphite under different conditions. Solid line: same as in the left panel; dotted 
line: supersonic approximation (Eq. |SJ; dashed line: neutral grain (or/and neutral gas); dotted-dashed line: (7=20- Volts charged grain 
moving in a T = 2 X 10^-K gas. 



sputtering is mostly due to the grain motion {non thermal 
sputtering) , as it can be seen in the right panel of Fig. |S| 
for V > 100 km s"\ the supersonic approximation (dotted 
line) is very close to the full calculation (solid line). In most 
cases, the grain potential acts to increase the effective sput- 
tering threshold, because part of the impact energy has to 
be spent in overcoming the Coulomb repulsion between the 
positively charged grain and the ions acting as projectiles. 
This is shown by the higher sputtering rate for a neutral 
grain for v ~ 100 km s"'^ (left panel: dashed line). If the gas 
temperature increases (left panel: dot-dashed line) the con- 
tribution of thermal sputtering becomes dominant for low 
velocities. 



3.4 Implementation 

The motion of dust grains is studied in a three-dimensional 
grid (typically made of 128^ cells). For each cell in the grid, 
gas density, temperature and ionization state are defined. 
After the determination of the position of the galaxies in the 
computational volume grid, grains are ejected into the IGM 
over the whole solid angle (from the center of each galaxy). 
The motion is assumed to start at a distance from the galaxy 
equal to its virial radius. A single ejected grain is thus char- 
acterised by its composition (graphite or olivine), radius a, 
velocity v and path direction. All these initial quantities are 
derived randomly with a Monte Carlo procedure from the 
adopted distribution (usually, those defined in § . Along 
the path the grain crosses cells with different gas conditions 
and its properties (velocity, charge and radius) change. 



As the grain enters a cell, it is instantaneously assigned 
an equilibrium potential (the time necessary to reach the 
equilibrium, t « \Z/Jce\ is usually much shorter than the 
time necessary for a grain to go across a cell; for the mean 
dust and gas conditions considered here, t ~25 yrs). For a 
fixed UVB, the equilibrium potential will depend on the gas 
density, temperature and plasma composition and on the 
grain material, radius and velocity. Given the equilibrium 
potential and the other gas and dust properties, we can com- 
pute the sputtering rate appropriate for the cell, dN/dt. The 
total number of atoms that are released as the grain moves 
through the ceU is then AN^{dN/dt){Al/v), where Al is 
the length of the grain path inside a cell. These sputtered 
atoms are deposited in the cell. As a result of sputtering, 
the grain loses mass (its radius reduces). 

Finally, we compute the reduction in velocity, Av = 
— (A//i')_Fdrag, due to coUisional and Coulomb drag. The 
grain velocity (and radius) is updated when the grain leaves 
a cell and enters the next along the path. The process is then 
repeated (computation of equilibrium charge, of deposited 
atoms and reduction in velocity) and the grain motion fol- 
lowed for the chosen time duration of the simulation t{ (un- 
less the velocity is so small that the grain will not move more 
than one cell in the remaining time). If Av is large in a cell 
the grain path is split in smaller units and Av, charge and 
AA'' are computed for each of them. We found that assuming 
lAwl/ii < 0.1 for each path segments provides accurate re- 
sults without overly increasing the computational time. For 
the cosmological simulation of path splitting is typically 
required at the end of a grain's path, for velocities below 100 
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km (the simulations of §0|are computed with a larger 
spatial resolution). 

In practice, the Monte Carlo procedure is run for A'^p 
cycles (or grain packets) to assure a good statistics (i.e. sam- 
phng of initial grain properties and path directions). If A^'g 
is the total number of grains that are ejected into the IGM 
in a simulation (defined by the mass of ejected dust and the 
size distribution), in each cycle A^c = Ng/Np identical grains 
will be ejected along a chosen direction. Thus, for each cell, 
Nc X AA'^ atoms are deposited in the IGM. The final out- 
put of the simulation is a grid with the number density of 
sputtered atoms. The final positions of the dust grains along 
their path is also stored. 



4 RESULTS: HOMOGENEOUS DENSITY 

To understand the details of grains motion and sputtering, 
we first analyze the case of a single galaxy ejecting grains 
into a homogeneous medium with z=3 mean cosmic den- 
sity. The IGM is assumed to be composed of H ll. He ll and 
electrons. Unless it is said otherwise, the gas temperature is 
assumed to be T = 2 x and the TRIM 2003 sputter- 

ing yield are used in the calculations. As we said in § 13.41 
we study the motion of the grain from the virial radius Ry, 
defined as the radius of the sphere e nclosing a mean den - 
sity of 200 times the critical density iNavarro et al.lll997D . 
In the cosmological simulation described in the next section, 
a galaxy of median mass has a dark matter halo of 1.5x 10^" 
M© (and a baryonic content of 1.5x10^ M0); the appro- 
priate value for the virial radius is Ry = 20 kpc (proper 
units). We follow the motion of each grain for ti=l Gyr: 
the results of the simulation can thus be compared to the 
observed properties of the universe at z = 2. During this 
time stretch, the UV B flux increases only by 25 per cent 
iBianchi et al 1 120011) : we neglect this time dependency in 
the calculations. 

It is instructive to study the velocity evolution of dust 
grains of different size and material. In Fig. |S]we show v as 
a function of the distance R from the center of the galaxy. 
Grains are ejected at Rv with three different initial velocities 
vo = 100, 500 and 1000 km sec~^. The two top panels refer 
to the case of grains of radius a=0.1-/im, while in the two 
bottom panels a—O.Ol-fim. Panels to the left are for graphite 
grains, those on the right for olivine. With the dashed line 
we show the results obtained for a grain whose charge is kept 
neutral (i.e. only the coUisional part is included in the drag 
force). For the same initial velocity, larger grains can attain 
larger distances from the galaxy: despite having larger cross 
sections when colliding with the IGM atoms, their decel- 
eration is smaller because they are heavier. For the same 
reason, heavier olivine grains travels to larger distances R. 
This dependence on the grain size and grain material den- 
sity can be derived analitycally for a neutral grain of radius 
a moving at supersonic speed. By integrating Eqn. |5| and 
|S|and neglecting the erosion of the grain due to sputtering 
(the radius a is kept constant), the velocity of the grain as 
a function of the distance R from the center of the galaxy 
can be written as 



The supersonic approximation for neutral grains of Eqn. 1111 
is shown as a dotted line in Fig. |^ (for ease of presenta- 
tion, only for iio=1000 km sec~^). The difference between 
the approximation and the curve for the case with only col- 
lisional drag is due to grain erosion. Despite small grains 
have smaller sputtering rates (see Eqn.|HJ they contain less 
atoms and their size reduces more. For uo = 1000 km sec~^, 
a=0.01-/im grains lose in the IGM approximately 3 x 10^ 
atoms, reducing their radius by 30-40 per cent (graphite 
and olivine, respectively). Conversely, a=0.1-/im grains lose 
2 X 10** atoms but their radius reduces less, by 20-25 per 
cent. 

The results for the complete calculation including grain 
charge and Coulomb drag is shown in Fig. |S| by the solid 
lines. When exposed to the z=3 UVB, 0.1-fim graphite 
grains attain a potential of about 25- Volts, while the poten- 
tial of olivine grains is about 20- Volts (the exact value de- 
pending on the grain velocity; 13. 2.111 . Because of Coulomb 
drag, the maximum distance R to which a grain can travel 
in the time t{ is reduced. Having a smaller charge (to- 
gether with the larger mass) , high velocity olivine grains can 
still travel at considerable distances (400-500 kpc) from the 
galaxy with velocities (kinetic energies) above the thresh- 
old for signiflcant sputtering. Instead, high velocity graphite 
grains eventually stop at R <300 kpc. Grains with a=0.01- 
^m attains potentials a couple of Volts larger than those 
with a=0.1-fj,m. As even in the neutral case small grains 
were not able to go beyond R «200 kpc, when the Coulomb 
drag is included they stop within 50 kpc from the ejection, 
even for the less restrictive case of a wo = 1000 km olivine 
grain. Small grains do not travel out to significant distances 
in the IGM. An analogous result has been obtained for the 
dust ejection from galaxies by some of the authors listed in 
§ 121 although with a more detailed study of the grain mo- 
tion, including gravitation and radiation pressure. This is 
why from now on we will use, as described in § |21 a flat size 
distribution, with 0.05< a[/im] <0.2, suggested by literature 
models. Together with this we will assume a fiat distribution 
of initial velocities, within the range 100< v <1000 km s~^. 

Fig-Qdisplays the carbon and silicon IGM metallicities^ 
resulting from the sputtering of carbon and silicon'^ atoms 
off graphite (left panel) and olivine grains (right panel) , for 
the adopted size and velocity distributions. The metallicity 
is reported as a function of the distance from the galaxy, R. 
The metallicity levels produced in our simulations depend 
linearly on the number of grains A'g ejected (equivalent to 
the ejected dust mass, for a fixed size distribution). The 
mass of dust ejected in Fig. |7| has been derived from the 
baryonic mass of our median z — 3 galaxy, by assuming 



v{R > R^ 



Vo exp 



3{R - Ry) 
Aap 



(11) 



^ The metallicity of element X is given by [X/H] = log(nx/nH) ~ 
log(nx/nH)0i where nx is the number density of element X 
and (X/H)q = log(nx/'^H)0 is the solar abundance. We use 
(C/H)q = -3.44 and (Si/H)© = -4.45 <Anders fc Grevessd 

^ As described in § 13.31 we have used the yield for the total num- 
ber of atoms sputtered off olivine. The number density of Si (and 
Fe, Mg) is thus 1/7 the total number density of atoms released by 
olivine grains in the IGM (the O number density being four times 
higher). The use of the rescaled total sputtering yield provides a 
good approximation to the specific sputtering yields for each of 
the species composing the material. 
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Figure 6. Velocity of dust grains injected into a homogeneous medium with mean 2 = 3 cosmic density, as a function of the distance 
from the source galaxy, R. For each panel, we show the velocity change of three grains with initial velocities (at Ry = 20 kpc) vq = 100, 
500 and 1000 km s~^ (from bottom to top). The two upper panels show the velocity of graphite (left) and olivine (right) grains of radius 
a=0.l-fim. The two lower panels show the same, but for a=0.01-/im. Solid lines show the results when the whole physics described in 
§|3]is taken into account, while dashed lines refer to the case in which the grain charge has been set to zero. For the vo = 1000 km s^^ 
we also show the solution for the neutral supersonic approximation of grains with a fixed radius fEgn. [TTl . The gas is made of Hll, Hell 
and electrons and has T=2 X lO^'-K. The TRIM 2003 sputtering yields have been used (SlOl. 



that, at maximum, 1/500 of the b aryonic mass is in dust in 
a galaxy jEdmunds fc Baled IToqI) and that 10 per cent of 
the total dust mass can reach the IGM (see Sect.j^J. In total, 
3x10* Mq of dust (both graphite and silicate) are ejected 
into the IGM for the simulation of Fig.Q Though we believe 
this to be a reasonable upper limit, the metallicity results 
can be easily scaled for any desired ejected mass of dust. 

We compare the results of our simulation w it h the 
met allicities of C and S i derived by ISchave et alj (l2003l) 
and lAguirre et alJ (l2004l) . They analyzed the pixel optical 
depth in a set of high resolution observations of the Lya 
forest using numerically simulated spectra and estimated 
the metallicity as a function of redshift and gas density. 
The gray areas in Fig. Q show the median metallicity (± 
la of lognormal scatter) for gas of mean density (overden- 
sity 5=1) at z = 2. The [Si/H] value has been extrapo- 



lated from the ratio [Si/C] jAguirre et alJl2004l) . with the 
caveat that it has been derived only for gas with 5 > 3. 
The scatter plotte d for [Si/H] is that derived for [C/H]. 
ISchave et alJ (|200^ find no significant trend for [C/H] ver- 
sus redshift in the range 1.8 < z < 4.1. A detailed study 
of the metal enrichment history is needed to asses whether 
this lack of evolution can be explained by dust sputtering or 
galactic winds (Aguirrc ct al. 2001c), or if instead an early 
enrichment from pr e-galactic {z ~ 9) sources is necessary 
jMadau et alJl2001^ . Here, we simply want to check if it is 
possible to reproduce, by using the dust sputtering mecha- 
nism only, the same level for [C/H] and [Si/H] as inferred 
from observations. 

Fig.Qshows that it is possible to obtain metallicity lev- 
els similar to those inferred from observations, at least for 
distances from the galaxy of R < 200 kpc. Within this range. 
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Figure 7. IGM metal pollution due to dust grains injected into a homogeneous medium with mean z = 3 cosmic density, as a function 
of the distance from the source galaxy, R. The left panels show the carbon metallicity [C/H] resulting from erosion of graphite grains, 
while the right panels show the silicon metallicity [Si/H] resulting from erosion of olivine grains. For each material, a flat distribution 
in the range 0.05< a [fJ,m] < 0.2 has been used for the grain radii and a flat distribution in the range 100 < r)o[km s~^] < 1000 for the 
grain initial velocities. The gas is made of Hll, Hen and electrons and has T=2 X 10'*-K. In total, 3 X 10^ Mq of dust has been ejected. 
The results shown with solid lines have been obtained using the TRIM 2003 sputtering yields, while for the dashed lines the Draine & 
Salpeter's yields have been used. The dotted lines refer to the TRIM 2003 case, but the gas temperature has been raise d to 2 X 10^-K. 
The gray areas show the metallicities derived at z = 2 for the overdensity 5 = 1 iSchave et al.ll2'o03l : lAguirre et ani2004) . 



the metallicity levels are quite similar for both materials. As 
we saw in Fig. |S] olivine grains can travel to larger distances 
and pollute with metals a larger area. However, the differ- 
ence in the trend can be appreciated only at _R > 300 kpc, 
where grains have reduced velocities and thus produce low 
levels of metallicity. We remind the reader that we have as- 
sumed, as an upper limit, the same relative abundance of 
graphite and silicate grains in the ejected dust. If silicates 
have a smaller ejection efficiency, as some studies suggest 
(see Sect.|5J the metallicity level in Fig. |7| will be lower (but 
the spatial distribution will be the same). 

A change in the gas temperature does not affect sig- 
nificantly the results. The dotted curve shows the case for 
r=2 X 10^-K. The metallicity becomes higher only for larger 
R, due to grains that are not supersonic anymore (see, e.g., 
the analogous case in Fig.|5J. Again, since graphite grains are 
slowed down more than those made of olivine, the difference 
is higher for [C/If] than for [Si/H]. Because the sputtering 
yield is quite flat for high energies, a further increase in the 
gas temperature does not produce a dramatic change in the 
results. Even for an extreme case in which T=2 x lO'^-K, the 
metallicities are only about 0.3-0.4 dex higher than those 
for r=2 X lO^'-K. In Fig. |7| we also show the results when 
the traditional Draine & Salpeter sputtering yields are used 
(dashed lines). The differences are small, especially when 
compared to the large uncertainties and scatter of the metal- 
licities inferred from observations. Graphite grains produce 
a lower [C/H] (but by only 0.2 dex) when using Draine & 
Salpeter yields, while the difference is smaller for olivine. 
This is mainly due by the difference in sputtering rates for 
V >200 km s~^, which is larger for graphite than for olivine 
(see the left panel of Fig. |KJ . For the rest of the paper, we 
will use the TRIM 2003 sputtering yields. 



5 RESULTS: COSMOLOGICAL SIMULATION 

We now study the motion of grains in a cosmological den- 
sity field. The simulation has been obtained using a multi- 
phase SPH code particularly designed for the study of galaxy 
formation (Marri & White 2003) and galactic winds. The 
initial parameters are those for the ACDM model adopted 
so far, {D.0 = 0.3, = 0.7, D.bh'^ = 0.028, as = 0.9). 
The simulation uses 128'' particles in a box of 10.5/i~^ co- 
moving Mpc. A group finding algorithm has been used to 
identify dark matter halos and the star-forming galaxies as- 
sociated to them. For a more d etailed descrip t ion o f the 
simulation, we refer the reader to lBruscoli et all (|200^ and 
iMaselh et ai] (l2004ll. 

We investigated the sputtering process on the z = 
3.27 simulation output. The particle properties have been 
mapped on a 3-D grid of 128^ cells. For each cell, we de- 
rived gas density, temperature and ionization fractions for 
each of the chemical species by performing a standard SPH 
smoothing on the 32 particles nearest to each cell. The spa- 
tial resolution of the grid is «82 kpc comoving (~27.5 
kpc physical). In total, 398 groups of particles have been 
identified as galaxies in the simulation box. The mass of 
dust in each galaxy has been derived from the baryonic 
mass, using the dust-to- baryonic mass ratio computed by 
lEdmunds fc Ealej Jl998ll . A gain, we allow 10 per cent of 
the dust to be ejected from the virial radius* As both the 



" In total, 2.7x10 Mq of dust are injected into the simulation 
volume, corresponding, for the adopted dust distributions, to a 
total number of dust grains Afg = 1.6 X 10^5. We have run sim- 
ulations with Np = 10^ packets of 1.6x10^^ grains each. This is 
sufficient for the results presented in this section to be indepen- 
dent of Np (i.e. Fig. ISlto ll2l are not affected by simulation cells 
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Figure 8. Cuts through the simulation box. The left panel shows the gas overdensity, while the centre and right panels show the carbon 
and silicon metallicities, respectively. The side of each map is 10.5h~^ comoving Mpc, 3.5 physical Mpc at z = 3.27. The unpolluted 
(white) cells in the [C/H] and [Si/H] maps are internal to the virial radius of the most massive galaxies in the selected slice. 



dust-to-baryonic mass ratio and the fraction of dust ejected 
are considered to be upper limit, we believe the metallicity 
results obtained in this paper (which scale with the number 
of ejected grains) to be upper limits to the dust contribution 
to the IGM pollution. The virial radius is computed from the 
mass of the halo associated to each galaxy, as specified in 
§ |1] Less massive galaxies are obviously more represented, 
with 80 per cent of the sample having a baryonic mass be- 
tween 8.5 X 10^ Mq (the minimum mass ) and 3.0 x 10^ M© 
(and a dark matter mass roughly 10 times that value). As 
we said in § |11 a galaxy of median mass ejects into the IGM 
3 X 10^ Mq of dust. Smaller objects also have virial radii 
which are quite close to the cell resolution (the halo of the 
smallest object has roughly the volume of a single cell). 

In Fig. IHl we show a cut through the whole simulation 
box, parallel to one of the box faces. The map in the left 
panel refers to the gas density field; the metallicity maps in 
the central and right panel ([C/H] and [Si/H], respectively) 
are the results of the calculations described in this paper. 
As discussed previously, the drag has a smaller effect on the 
heavier and less charged silicate grains. This is evident in 
the larger area of the slice polluted with silicon. The slice 
passes through the most massive galaxy in the simulation 
(with baryonic mass 7 x 10" Mq) which can be identified 
in the metallicity maps with the largest "hole" (inside a 
high metallicity region). The holes are the regions within 
the halo of each galaxy, that are excluded from the simula- 
tion (as grains are assumed to move from the virial radius 
outward). For the largest galaxy, i?v=70 kpc physical (2.5 
times the dimension of a cell). Other high mass objects can 
be seen to the right. It is interesting to note how the massive 
galaxies are not necessarily surrounded by a wider metal en- 
riched area, despite injecting into the IGM a larger amount 
of dust. As they reside in denser gas, grains ejected from 
these objects are more easily stopped by the drag. 

In Fig. 1^1 we show the metallicity as a function of the 
gas overdensity. Each dot in the left (right) panel represents 
a simulation cell which has been polluted by carbon (silicon) 

where only a few grain has passed and thus not having converged 
to a mean value for the metallicity, grain density and number). 



as the result of the passage of one or more graphite (silicate) 
grains. As a reference, in the left panel we have plotted the 
median carbon metall i city ( ± la of lognormal scatter) de- 
rived bv lSchave et al.l i2003l) as a function of the gas over- 
density (gray area). We are showing the value at z = 2. In 
the right panel, the same are a is indicated, but sc aled for 
the [Si/C] ratio measured bv lAguirre et all (120041) for gas 
with 5 > 3. S i milar metallicity levels have been reported by 
ISimcoe et al.l J^QO^) , by fitting hydrogen and metal lines at 
2 = 2 — 2.5 in QSO's absorption spectra. They do not de- 
tect a significant trend with the gas density. However, this 
may be due to the different choice for the UV background 
needed to correct for ioniza tion, as harder spec tra lead to a 
shallower dependence on 6 JSchave et alJl2003h . 

As seen in the homogeneous case, grain sputtering 
can produce non negligible metallicities in the IGM gas. 
The dependence of metallicity on density is steeper than 
what inferred from observations. Most cells with 5 « 1 
have [C, Si/H] « -5, approximately one order of magni - 
tud e lower than th e value s measured bv lSchave' et al.1 ll2003h 
and lAguirre et al.l ll2004l) (especially [C/H]). The contribu- 
tion of dust sputtering becomes more important for mod- 
erately overdense gas with 5 = 10 — 100, which can be 
translated in the neutral hydrogen column density range 
14.5 < logA''(Hl) < 16. We find, however, that only 4-5 
per cent of the cells in this overdensity range have [Si, C/H] 
values within the observationally allowed area. We remind 
that we believe this to be an upper limit, both for the total 
amount of dust and for the relative number of silicate and 
graphite grains ejected into the IGM (we use a 1:1 ratio). 
If the fraction of silicate grains is smaller, as some works 
seem to suggest (§ [Si/H] would be smaller. However, 
the trend with density will be the same. While the level of 
metallicity in our simulations scales with the ejec ted mass 
of dust, the value derived bv ISchave et alJ i2003l) depends 
on their assumption for the UV background. The scatter, 
instead, depends less on the assumptions and more on the 
underlying physics and on the density structure of the IGM 
(Aguirre, private communication). For cells with 5 « 1, we 
find a log- normal scatter of about 1.6 dex (slightly larger for 
Si). This is about twice what inferred from observations. 
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Figure 9. Mctallicity vs overdensity for simulation cells polluted by carbon from graphite grains (left) and silicon from olivine grains 
(right). Each cell is indicated by a dot. Contours are plotted for regions where the density of dots becomes high. In both panels, contours 
start at 20 (and inc reas e by 40) cells per 0.1 dex in overdensity and metallicity. The gray area refers to the metallicity measurements of 
ISchave et all hOO^ and lAguirre et al.l i2004) . 



The features in the scatterplot reflect the movement 
of dust grains. In each panel, the region in the top-right is 
occupied by cells just outside the virial radius of the injecting 
galaxies, where the grain movement st arts. In principle, fo r 
an isolated halo, gas at _Rv has S ~60 iNavarro et al.lll997t) . 
However, due to the complex IGM structure and to the lack 
of resolution, the gas cells where the grain movement starts 
have overdensities 1 < log]^o(^) 3 (we will comment on 
resolution in § (SJl. Because of the high density and grain 
velocities, metallicities (and sputtering rates) are high in 
these cells. As the grain moves towards the less dense gas, it 
is slowed down by gas drag. For reduced velocities and gas 
densities, the sputtering rates are lower and so the resulting 
metallicities. This explains the general trend from high p - 
high metallicity to low p - low metallicity seen in both panels 
of Fig. 1^ The slope of the relation depends both on the 
initial grain velocity and radius (as seen for the homogeneous 
case of § 0]l and on the gas density. 

There are some exceptions to such trend. Grains ejected 
into the higher density gas are stopped by drag more effi- 
ciently and pollute a smaller region around them. This is 
why the high metallicity cells at 5 > 2 appear isolated (i.e. 
with no contiguous lower metallicity and S points) in the 
scatterplot. As the velocity quickly reduces within the dense 
region, some high-(5 cells end up with low metal pollution. 
In this case the slope of the relation becomes very steep 
and lack of resolution limits the detailed study. Other low 
metallicity cells with high-5 mark the end-point of grains 
previously travelling in a less dense environment that stop 
as they cross a denser filament in the density field. 

The gas which grains encounter as they leave a galaxy 
is often shock-heated by gravitational collapse and/or SNe. 
The cells with the highest metallicities in Fig.|n|have temper- 
atures in excess of 10^-K. However, this has little effect on 
the metallicity levels: the grain initial velocities considered 
here are always higher or close to those delimiting the non- 
thermal regime for sputtering. Instead, the high temperature 



determines the mechanism responsible for the drag. In the 
hot gas, collisional drag is more important 13.111 and the 
grain movement does not depend on the charge. Coulomb 
drag becomes more important than collisional drag in the 
low temperature gas close to 5=1 and it is because of it that 
grains come to a halt. Within the adopted simulation time 
tt, most grains are slowed down and have energies below the 
threshold energy for sputtering. Even a ten-fold increase in 
tf does not change the results significantly. 

By counting cells which have been polluted, we can pro- 
vide an estimate for the volume filling factors of metals. A 
wider volume is polluted by silicon, and this is clearly shown 
by the larger number of points (cells) in the right panel of 
Fig. |U] (and by the map in Fig. |5J . By checking simulations 
with increasing number of grain packets, we found that the 
silicon filling factor tends asymptotically to about 18 per 
cent. As most of the cells are occupied by low density gas far 
from sources, filling factors are obviously higher for denser 
regions. When considering all cells with overdensity S > 10, 
the silicon filling factor increases to 40 per cent. The mean 
filling factor of cells polluted with carbon is a factor three 
lower (~6 per cent). Thus, the distribution of metals in our 
simulations is rather inhomogeneous. Though a fraction of 
cells in Fig.Elhave metallicit ies compatible with the median 
value of ISchave et alJ i2003l) . the median metallicity in our 
simulation is zero, for any gas density bin. 

As Fig.|^clearly show, the ratio [Si/C] in our simulation 
is lowe r than 0.77±0.05, the value inferred by^guirrc ct ^ 
ll2004h at 2; « 3 for gas with 3 > 3. In the high density 
gas, [Si/C] depends mostly on the assumption about the 
relative number of silicate/graphite grains and on the dif- 
ferences between the sputtering yields. If a cell is crossed by 
the same amount of both kind of grains, and if the kinetic 
energy is above the sputtering threshold, the ratio is quite 
constant. For example, silicate and graphite grains travel- 
ling at 500 km s^^ in the same cell would pollute the gas to 
[Si/C]~0.15. Indeed, the high density cells where the grain 
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Figure 10. Initial (solid line) and final (dashed line) grain size 
distributions. The two histograms have been normalised to unity. 
The distribution for grains at rest in gas with 5 < 10 is shown 
with a dotted line (about 10 per cent of the total). 

movement starts have [Si/C]~ 0.1 ± 0.08. Again, a com- 
parison between the scatter in simulations and that derived 
from observations may be less dependent on the assump- 
tions that have been made: for the high density cells, the 
scatter is similar to that of lAguirre et alj i2004l) . while it 
increases at lower density. The [Si/C] value increases as well 
at lower S, because olivine grains keep higher velocities to 
larger distances (for our assumptions, 3 per cent of the sim- 
ulated volume has [Si/C]> 1). 

During their travel, dust grains are eroded but never 
destroyed by sputtering. At the end of the simulation, dust 
grains at rest are present throughout all the volume pol- 
luted with metals (the filling factors of metal polluted cells 
and cells containing grains at rest are similar). The final 
size distribution (Fig. llOt is roughly flat, as the input one, 
but with a different range in radii. The smallest grains have 
radii a ~0.02-/im and originate from particles with initial 
radius a—0.05-^m, while largest grains reduce their radius 
from a=0.2-/im to a=0.15-/im^. The latter are the grains 
that contribute mostly to metal pollution and, as we dis- 
cussed previously, travel farther out from their injection 
point: for regions of diffuse gas {S < 10) the grain distri- 
bution is peaked around sizes 0.1-0.15 /xm. There is no sig- 
nificant difference between graphite and silicate grains. 

In Fig. 111! we show the final grain density (both silicate 
and graphite) as a function of the gas overdensity, for cells 
occupied by dust grains at rest. Features corresponding to 
those of Fig. ^ can be noticed. The higher grain densities, 
corresponding to higher gas densities, are due to a relatively 
large number of grains that stop close to the injection point. 
Only a limited number of grains travel far out, and this 
ma kes for the reduction of g rain density for 5 < 10. 

Ilnoue fc Kamaval (l2004l) derived an upper limit for the 
dust grain density. After adopting various cosmic star for- 
mation histories (to account for the dust generation) they 

9 Incidentally, by using TRIM iGrav fc Edmundj |20^ have 
found that sputtering does not simply lead to grain erosion: the 
stopping of projectiles into the material may change the structure 
and composition of the grain, especially in a metal rich gas 
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Figure 11. Dust grain density in the IGM vs gas overdensity for 
simulation cells occupied by grains at rest. Each cell is indicated 
by a dot. Contours are plotted for regions where the density of 
dots becomes high; they start at 20 (and increase by 80) cells per 
0.1 dex in over density and grain density . The upper limit refers 
to the results of llnoue fc Kamaval 1200^^) . 

constrained the amount of dust in the IGM at 5 = 1 with the 
maximum extinction and reddening allowed by the (now pre- 
ferred) cosmological explanation for the dimming of high-z 
Type la SNe. Furthermore, they tightened the constraint re- 
quiring that dust grains cannot significantly alter, via photo- 
electric heating, the thermal history of the IGM derived from 
the l ow column density lines of the Lya forest iSchave et al] 
|2p00). The upper limit for grains of size a = 0.1-/xm is plot 
in Fig. 111! The median density for 5 = 1 cells occupied by 
grains is about 1.5 orders of magnitude lower than the upper 
limit. As for metallicity, the results of Fig. Illl scale linearly 
with the amount of dust ejected in the IGM. If our upper 
limits are correct, the extinction and heating effects of IGM 
grains are not likely to be detected. 



6 COMPARISON WITH PREVIOUS WORKS 

lAguirre et alj i2001al l9) found that dust expulsion and ero- 
sion can account for the observed levels of carbon and sil- 
icon enrichment in the IGM. Since there are several differ- 
ences between their work and ours, it is not easy to com- 
pare the results. They position grains at a distance from the 
galaxy where gravitation is balanced by radiation pressure. 
Because of this, silicate grains (having a larger material den- 
sity) do not lay as far out as graphite grains. Furthermore, 
silicate grains have a smaller efficiency for radiation pres- 
sure. This is opposite to our results, where heavier silicate 
grains can travel to larger distances. They adopt essentially 
a power law distribution for grain sizes. Grains with radii 
a = 0.03 — 0.05/im have larger equilibrium distances than 
grains with a > 0.1-^m. Again, this is opposite to what we 
find here by studying the dynamics of grains. 

An other differenc e is in the amount of dust that is 
ejected. lAguirre et alJ allow half of the metal content of each 
galaxy in their simulation to be distributed in the IGM as 
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dust. As typically only half of a galaxy's metal mass is be- 
lieved to be in dust, this is equivalent to ejecting the whole 
dust content from the galaxy. Although we believe that 10 
per cent is an upper limit to the fraction of dust ejected 
into the IGM (see even allowing all the dust mass to 

travel to the IGM (with large grains overrepresented with 
respect to the standard MRN power law) does not change 
significantly the results. As it can be seen by simply shift- 
ing upward by 1 dex the metallicities of Fig. El more metal 
polluted cells will lay within the range derived from observa- 
tions. Still, the filling factor of metal polluted cells will not 
change. The same shift applies to the points of Fig. 1111 Even 
increasing by one order of magnitud e the grain density, the 
limit set bv llnoue fc Kamaval 11200^ will not be violated in 
most cells with grains at rest. 

The distribution of metals in our simulations is inho- 
mogeneous, with filling factors smaller than 0.5 for any gas 
density bin, and it seems smaller than that in lAguirre et alJ 
i2001al l3) (we compare our results to their case in which 
only the metals effectively eroded from grains are consid- 
ered). One potential interpretation is that in our simulation 
the IGM density field has been "frozen" at 2 = 3, where all 
grains have been ejected. The 40 per cent decrease of the 
mean density from z = 3 to z = 2 (the end of the simula- 
tion) it is however unlikely to produce significantly different 
effects. The lowest mass galaxies in our simulations have ha- 
los of the same size as the cells. Lacking higher spatial (and 
density) resolution, grains ejected from small galaxies at the 
virial radius may find themselves in a denser environment 
than expected. As both these issues concern the movement 
of grains in less dense environments, we ran a simulation in 
which grains were injected at a distance from each galaxy 
twice the virial radius. The metallicity-density plot for this 
test simulation is shown in Fig. 1121 Now the cells where the 
grain movement starts have a slightly lower density (there 
are less atoms to collide with) and obviously have lower 
metallicity (the peak of its distribution is shifted towards 
S ~1). The reduction in gas particles also mean a smaller 
drag and a smaller reduction in velocity as the grain moves. 
The slope of the metallicity-densi ty trend in Fig . lT^ becomes 
flatter, and closer to the data of ISchave et al.l (2003). Also, 
a good fraction of J « 1 cells has now metallicities within 
the shaded area (at least for [C/H]). There is also an obvious 
increase in the filling factors, with silicon being present in 50 
per cent of the volume, and carbon in 20 per cent. Similar 
trends can be seen in the grain density distribution, with 
an increase of the density by about 0.2 dex for cells with 
5 < 1 (a reduction by 0.2 dex for 5 > 10). Clearly, a better 
spatial resolution is desirable, although the conclusions are 
not going to be substantially altered. 

Finally, we note here that it would be possible to obtain 
higher metallicity levels by destroying completely the grains, 
via some unknown mechanism other from thermal/non- 
thermal sputtering. In our simulation a sizable fraction of 
large grains are at rest in gas with 5 < 10. By adding their 
contribution in atoms to the metallicity produced by the 
sputtering processes we have considered, we would rise by 
almost 2 dex the metallicity of one third of cells at 5 < 1, 
while leaving metallicity levels unaltered for S > 10. 



7 SUMMARY AND CONCLUSIONS 

We have analysed the motion of dust grains in the z — 3 
IGM, to study the erosion due to sputtering and the subse- 
quent pollution of the gas with metals. A fraction of the dust 
grains in a galaxy are thought to escape the gravitational 
well because of radiation pressure from starlight. Using re- 
sults in the literature, we have defined plausible distribu- 
tions for the sizes and velocities of the escaping grains. Our 
fiducial grain size distribution is fiat, with radii in the range 
0.05 < a[/im]< 0.2, as literature results suggest that only 
large grains can escape from galaxies. The assumed velocity 
distribution is flat as well, with 100 < v[km s~^]< 1000. 
For such velocities, grains are supersonic under most of the 
gas conditions explored: thus, non-thermal sputtering is the 
dominant mechanism for the grain erosion. 

The motion of each ejected grain is studied from the 
virial radius, until the velocity becomes small and the grain 
stops. Our computation takes into account: 

(i) a calculation of the grain charge, due to collisions with 
electrons and ions in the IGM plasma and to photoejection 
of e lectrons by a metaga lactic UV background. We have used 
the lSianchi et afl 1I2OO j) UV background including the con- 
tribution of galaxies and taken into account the velocity de- 
pendence of the coUisional charging rates; 

(ii) the gas drag, including both the coUisional term and 
the Coulomb term due to interaction of the (heavily) charged 
grain in the z — 3 gas and the charged particles in the gas; 

(iii) the sputtering of atoms off the grain surface, due to 
collisions with the ions in the IGM. Using the code TRIM, 
we have derived new sputtering yields for H and He atoms 
colliding with graphite and a silicate, olivine. Analytical fits 
to these yields are provided. The new yield have been com - 
pared to the widely used yields of lDraine fc Salpeterl il979t) . 

The grain motion has been studied for a single source 
of grains (galaxy) in a homogeneous density field and inside 
a cosmological simulation, allowing the galaxies detected in 
the simulation box to eject 10 per cent of their dust mass 
into the IGM. We believe this assumption (as well as the re- 
sulting metallicities and grain density) to be an upper limit, 
compatible with the presumable MRN-like dust distribution 
inside a galaxy and the adopted size distribution for ejected 
grains. The motion is followed for 1 Gyr (half the Hubble 
time at z = 3), allowing us to compare our results to obser- 
vations at 2 = 2. The main results are: 

(1) Only large (a > 0.1-fim) grains can reach out for the 
low density regions in the IGM, because the drag is less ef- 
fective for heavier grains. For the same reason, silicate grains 
are able to diffuse in the IGM more than graphite grains and 
pollute with metals a larger fraction of the cosmic volume. 
Carbon and silicon metallicities produced by dust sputtering 
of graphite and olivine grains show a well defined trend with 
gas density. High density cells in the vicinity of galaxy halos 
(where the grains start their travel) have high metallicities, 
whereas fewer grains (originally at high v) reach gas with 
5 ~1, with reduced velocity producing smaller sputtering 
rates. Those grains that are able to reach lower density re- 
gions are stopped, within the simulation time, by Coulomb 
drag. None of the large grains is completely destroyed by 
sputtering. The results summarised in this point are general 
and do not depend on the assumption on the adopted size 
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and velocity distribution and fraction of dust mass ejected 
from galaxies, nor on the specific density field used in this 
work. 

(2) For our parameter choice, most of the 5 ~ 1 simu- 
lation cells polluted with metals have [C, Si/H] ~ -5, ap- 
proximately one order of magni tude lower than the me- 
dian values r e cently derived by ISchave et al.l ll2003ll and 
lAguirre et alj i2004) . The contribution of dust sputtering 
becomes more important for moderately overdense gas with 
S = 10 — 100, corresponding to neutral hydrogen column 
densities 14.5 < log N{¥^ l) < 16. For this overdensity range, 
however, only 4 per cent of the cells have [C/H] values within 
the observationally allowed area. For any density bin, metal- 
licities have a large scatter (a lognormal a > 0.15 dex at 

5 Ril). This is twice what has been inferred from observa- 
tions of absorption lines. In addition, we recall that olivine 
grains pollute the IGM with Fe, Mg (same level as [Si/H]), 
and oxygen ([O/H] is a factor four higher than [Si/H]), al- 
though these species are less readily observed in the IGM. 

(3) Polluted cells in the range where dust sputtering is 
important (10 < 5 < 100) tend to have 0<[Si/C]<0.2, al- 
most independently of the gas density. The value we have 
obtained strongly depend on the graphite/silicate propor- 
tion, which we have assumed to be 1:1, as for MRN dust. If 
this is correct, we would have [Si/C]<0.2 when allowing for 
a lower ejection efficiency of silicate grains with respect to 
graphite (as suggested by a few work in the literature) . Such 
ratio is substantially smaller than the nucleosynthetic one, 
approximately equal to 0.9. Hence, this ratio might repre- 
sent a valuable tool to assess the importance of the grain 
sputtering enrichment by systematic searches and studies of 
low [Si/C] absorption systems. 

(4) Resulting metal distributions are very inhomoge- 
neous, with only 18 per cent of the volume occupied by cells 
polluted with silicon (from far-reaching olivine grains) and 

6 per cent with carbon. These results slightly depend on 
resolution, but we estimate that filling factors larger than 
50 per cent (for silicon) for cells at 5 = 1 enriched by dust 
sputtering are unlikely. 



(5) Even allowing all dust mass produced inside galax- 
ies by 2 = 3 to be ejected, the grain density in the low 
(S ~ 1) density IGM i s below the upper limit derived by 
llnoue &i Kamaval ll2004f) , from IGM thermal history and ex- 
tinction considerations. However, some moderately dense re- 
gions with (5 > 10 might contain a significant amount of dust, 
whose extinction effects remain to be calculated in detail. 

To conclude, our results show that radiation-pressure 
driven dusty fiows, followed by sputtering of the grains, 
might represent a viable and attractive mechanism to enrich 
the IGM with (cool) heavy elements. This process mostly af- 
fects moderately overdense IGM, without affecting too much 
the low-density gas: this is understood by the enhanced sput- 
tering efficiency where the gas density is higher. 

At least, three aspects of the problem need further study 
before we can draw final conclusions. The first concerns the 
ejection of dust grains from galaxies: although a few au- 
thors have studied the problem (we have used their results 
here as initial conditions) so far no work presents a detailed 
statistical analysis of the properties of escaping grains. In 
particular, estimates are lacking for the fraction of the dust 
mass which is ejected, as a function of the dust size and 
composition and of the galactic mass. This is important if 
we want to study the history of dust ejection and metal pol- 
lution by grain sputtering. At reddhift z > 3, despite an 
higher IGM density (increase in the drag), the physical dis- 
tance between galaxy is smaller and grains may be able to 
pollute a larger volume of the universe. However, galaxies 
have a smaller stellar mass: if the internal radiation field is 
weaker, the dust ejection efficiency may be reduced (unless 
this effect is contrasted by changes in the ISM drag and in 
gravity). The second aspect relates to the UV background 
fluctuations: in this study we have considered the UVB to be 
isotropic. However, radiative transfer effects (as shadowing 
and shielding), in the vicinity of galaxies or dense cosmic 
filaments, might create anisotropies resulting in a net radi- 
ation force on the grain. Such force might be large enough 
to prevent the grain from stopping as it currently happens 
in our simulations. In this case, the grain can continue to 
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move essentially taking a random walk, until it enters a re- 
gions of very high density (for example an accreting flow) 
within which it becomes fully coupled to the gas. If so, more 
metals could be released and less grains would survive, thus 
further decreasing the resulting intorgalactic extinction. The 
third important effect concerns the intergalactic magnetic 
fields, which we have so far ignored. Although they might 
not be generally dynamically important, they might affect 
dust sputtering, via the betatron effect. We plan to focus on 
these additional physical effects in forthcoming papers. 
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